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Abstract. 

We study the accuracy with which the WIMP mass could be determined by a 
superCDMS-like direct detection experiment, given optimistic assumptions about the 
detector set-up and WIMP properties. We consider WIMPs with an interaction cross- 
section of cr p = 10~ 7 pb (just below current exclusion limits) and assume, initially, that 
the local WIMP velocity distribution and density are known and that the experiment 
has negligible background. For light WIMPs (mass significantly less than that of the 
target nuclei) small variations in the WIMP mass lead to significant changes in the 
energy spectrum. Conversely for heavy WIMPs the energy spectrum depends only 
weakly on the WIMP mass. Consequently it will be far easier to measure the WIMP 
mass if it is light than if it is heavy. With exposures of £ = 3 x 10 3 , 3 x 10 4 and 
3 x 10 5 kg day (corresponding, roughly, to the three proposed phases of SuperCDMS) 
it will be possible, given the optimistic assumptions mentioned above, to measure the 
mass of a light WIMP with an accuracy of roughly 25%, 15% and 2.5% respectively. 
These numbers increase with increasing WIMP mass, and for heavy WIMPs, m x > 
0(500 GeV), even with a large exposure it will only be possible to place a lower limit 
on the mass. Finally we discuss the validity of the various assumptions made, and 
the consequences if these assumptions are not valid. In particular if the local WIMP 
distribution is composed of a number of discrete streams it will not be possible to 
determine the WIMP mass. 
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1. Introduction 

Diverse cosmological observations indicate that the majority of the matter in the 
Universe is dark and non-baryonic (e.g. Ref. [I]). Weakly Interacting Massive Particles 
(WIMPs) are one of the leading cold dark matter candidates, and supersymmetry 
provides a concrete, well-motivated, WIMP candidate in the form of the lightest 
neutralino (e.g. Ref. [2j [3]). The direct detection of WIMPs in the lab [1] would not 
only directly confirm the existence of dark matter but would also probe the parameters 
of supersymmetry models. Constraints on, or measurements of, the WIMP mass and 
interaction cross-section will be complementary to the information derived from collider 
experiments [5], [6] . It is therefore pertinent to examine the accuracy with which WIMP 
direct detection experiments will be able to measure the WIMP mass, if they detect 
WIMPs. 

Direct detection experiments could potentially measure the WIMP mass via the 
mass dependence of either the energy spectrum [7j or the 'crossing energy' at which 
the phase of the annual modulation signal [8j, due to the motion of the Earth, changes 
sign [9], [lOl [H]. The size of the annual modulation signal is small (of order a few per- 
cent), however, and given the current exclusion limits from CDMS [T2l [T3] it is unlikely 
that even planned tonne scale experiments, such as SuperCDMS [HI US], Xenon [16] 
and EURECA [17J, will be able to accurately measure the energy dependence of the 
annual modulation phase [101 [18] 0. Furthermore the shape, phase and amplitude of 
the annual modulation signal depend sensitively on the detailed local WIMP velocity 
distribution [22] whereas the shape of the differential event rate does not [23j [21], if 
the local WIMP distribution is smooth. The mass dependence of the differential event 
rate [7] (see also [25]) therefore appears to offer the best prospect, in the short to medium 
term at least, for probing the WIMP mass using direct detection experiments. 

In this paper we examine the accuracy with which a future SuperCDMS [El [15] like 
direct detection experiment will be able to measure the WIMP mass, given a positive 
detection. We consider the effect of varying the underlying WIMP mass and detector 
exposure, and also examine the uncertainties which arise from our lack of knowledge of 
the underlying WIMP distribution. In Sec. [2]we outline the calculation of the differential 
event rate and its dependence on the WIMP mass. We describe our Monte Carlo 
simulations and results in Sec. [3], discuss the validity of the assumptions made in Sec. 0] 
and conclude with discussion in Sec. [51 

| The DAMA collaboration have, with a Nal detector and an exposure of ~ 1.1 x 10 5 kg day, observed 
an annual modulation, which they interpret as a WIMP signal 19 , however it appears to be possible to 
reconcile this with the CDMS exclusion limit only by invoking 'non-standard' WIMP properties (such 
as the WIMP-proton and WIMP-neutron couplings being different [5D] or WIMPs which can scatter 
inelastically [21]) . 
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2. Differential event rate 

2.1. Basic calculation 

The direct detection differential event rate, or energy spectrum, depends on the WIMP 
mass, its interaction with the detector nuclei and the velocity distribution of the 
incoming WIMPs. Assuming purely spin-independent coupling, the event rate per unit 
energy, is given by (see e.g. [2j [7j): 

^^^(W, (i) 

where p x is the local WIMP density, cr p the WIMP scattering cross section on the 
proton, yU PX = {m v m x )/{m p + m x ) the WIMP-proton reduced mass, A and F(E) the 
mass number and form factor of the target nuclei respectively and E is the recoil energy 
of the detector nucleus. The dependence on the WIMP velocity distribution is encoded 
in J-"(E), which is defined as 

HE) = ( ^^dv) , (2) 



where f E (v,t) is the (time dependent) WIMP speed distribution in the rest frame of 
the detector, normalized to unity and (..) denotes time averaging. This is calculated 
from the velocity distribution in the rest frame of the Galaxy, / G (v), via Galilean 
transformation: v — > v = v + v E (t) where v E (t) is the Earth's velocity with respect to 
the Galactic rest frame. The lower limit of the integral, t> m m, is the minimum WIMP 
speed that can cause a recoil of energy E: 

1/2 

(3) 



Em A \ 



where uia is the atomic mass of the detector nuclei and [ia x the WIMP-nucleon reduced 
mass. We use the Helm form factor [26] with parameter values as advocated by Lewin 
and Smith [TJ. For most calculations we use the 'standard halo model', an isotropic 
isothermal sphere, for which the local WIMP velocity distribution, in the Galactic rest 
frame, is Maxwellian (c.f. Ref. [8]) 

/G(v) = N [exp (-|v| 2 /^c) - exp (-r&>c)] l v l < , (4) 
/ G (v) = |v|>*w, (5) 

where iV is a normalization factor, v c = 220 km s -1 |27J and v esc = 540 km s" 1 [28J are 
the local circular and escape speeds respectively and we use the usual fiducial value 
for the local WIMP density, p x = 0.3GeVcm -3 . We use the expressions for the time 
dependence of the Earth's velocity with respect to the Galactic rest frame from Ref. [7], 
to transform the velocity distribution into the lab frame. 
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Figure 1. The differential event rate in the E — > OkeV limit, (dR/dE)o, (top 
panel) and the characteristic energy scale, E'r, (bottom panel) as a function of the 
WIMP mass, m x for a Ge detector. The solid (dotted) lines are the fits to the full 
calculation including the Earth's velocity and the Galactic escape speed for threshold 
energy E'th = (10) keV. The short dashed lines are eqs. (0 and ([8]) which neglect the 
Earth's velocity and the Galactic escape speed. In the bottom panel the long dashed 
lines are the asymptotic mass dependences of Er in the m x <C wia and 3> wa limits 
(£r oc m x and cx const respectively). 



2.2. Dependence on the WIMP mass 

The WIMP mass dependence of the differential event rate can most easily be seen, 
following Lewin and Smith |7J, by first neglecting the Earth's velocity and the Galactic 
escape speed. In this case eq. ([I]) can be written as 

™ {E) =m\^(e\ nE) , (6) 



dE y ' \dEj r V E R/ 
where (dR/dE) , the event rate in the E — > OkeV limit, and the characteristic 
energy scale, are given by 



dEJ ^nl x m x v c 
and 

^H = ^, (8) 

respectively. When the Earth's velocity and the Galactic escape speed are taken into 
account eq. flTJ is still a reasonable approximation to the event rate, provided that 
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multiplicative constants c and cg R , are included in the expressions for (^§) and £r [7j. 
The exact values of the constants, which are of order unity, depend on the target nuclei, 
the energy threshold and the Galactic escape speed. We find, using least squares fitting 
to the full numerically calculated time averaged differential event for a Ge detector with 
energy threshold E t h = OkeV, Co ~ 0.78 and ce b ~ 1.72, with a weak dependence (in 
the 3rd significant figure) on the WIMP mass. 

The mass dependence of (dR/dE) and E R for a Ge detector are shown in fig. [TJ 
As expected from eq. (171), (dR/dE) oc m" 1 , since m x 3> m p . Similarly as expected 
from eq. ([8]), for m x < m^, E-& oc m 2 . This mass dependence weakens with increasing 
WIMP mass, and Er tends to its constant large mass asymptote for m x ~ 0(1 TeV). 
We will see below that this has important consequences for the determination of the 
WIMP mass from direct detection experiments. If, anticipating section III, we only 
fit to the event rate above an energy threshold E t h = lOkeV we find that for WIMP 
masses below ~ 100 GeV the fitting constants depend on the WIMP mass (this is visible 
in fig. CD as the deviation of the dotted lines from the solid lines for small WIMP masses). 
This illustrates that while eq. [6] is a useful approximation for demonstrating the mass 
dependence of the energy spectrum, for concrete applications the time average of the 
local WIMP velocity distribution should be calculated explicitly. 

3. Monte Carlo simulations 

We use Monte Carlo simulations to examine, for a range of detector exposures and 
input WIMP masses, how well the WIMP mass could be determined from the energies 
of observed WIMP nuclear recoil events. 

3.1. Detector properties 

Our simulated detector is based on the proposed SuperCDMS experiment [THUS], being 
composed of Ge with a nuclear recoil energy threshold E t ^ = 10 keV. We assume that 
the background event rate is negligible, as is expected for this experiment located at 
SNOLab [Tl"l[T5] . and that the energy resolution is perfect. For simplicity we assume that 
the nuclear recoil detection efficiency is independent of energy. The energy dependence 
of the efficiency of the current CDMS experiment is relatively small (it increases from 
~ 0.46 at E = 20keV, the energy above which data from all detectors is analyzed, to 
~ 0.51 at 100 keV pjlj). For further discussion of these assumptions see Sec. EH 

The proposed SuperCDMS detector consists of 3 phases, with detector masses of ~ 
25 kg, 150 kg and 1 ton. We consider efficiency weighted exposures [§| S = 3 x 10 2 , 3 x 10 3 , 
3 x 10 4 and 3 x 10 5 kgday. The later three exposures correspond, roughly, to a detector 
with mass equal to that of the 3 proposed phases of SuperCDMS taking data for a year[JJ] 
with a ~ 50% detection efficiency. 

§ For brevity we subsequently refer to this as simply the exposure. 

|| Accumulating this much 'live-time' would of course take substantially longer than a year. 
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3.2. WIMP properties 

We assume a fixed WIMP-nucleon cross-section of er p = 10~ 7 pb, which is just below 
the current exclusion limit from the CDMS experiment [12]. Since the number of 
events expected in a given experimental set-up is directly proportional to the cross- 
section, there is a straight-forward scaling to other cross-sections i.e. an exposure 
£ — 3 x 10 4 kgday for a p = lCT 7 pb is equivalent to an exposure of £ = 3 x 10 5 kgday 
for cr p = 10~ 8 pb. We consider input WIMP masses of m x = 25,50,100,250 and 
500 GeV, with a Maxwellian speed distribution with circular speed v c = 220 km s -1 . For 
m x = 100 GeV we also consider circular speeds in the range v c = 180 to 260km s -1 [2Z] 
and variations in the form of the speed distribution. 



3.3. Statistical analysis 

We estimate the WIMP mass and cross-section by maximizing the extended maximum 
function (which takes into account the fact that the number of events observed in a 
given experiment is not fixed) e.g. Ref. |29j : 

A^p'exp (-A). 



-U^rf(Ei) • (9) 



Here N exp t is the number of events observed, E- X (i = 1, N expt ) are the energies 
of the events observed, f(E) is the, normalized, differential event rate and A = 
£ {dR/dE) dE is the mean number of events (f(E) and A depend on m x and er p ). 
We calculate the probability distribution of the maximum likelihood estimators, for each 
exposure and input WIMP mass, by simulating 10 4 experiments. We first calculate the 
expected number of events, A in = £ JJ° (dR/dE) dE, from the input energy spectrum. 
The actual number of events for a given experiment, N expt , is drawn from a Poisson 
distribution with mean A in . We Monte Carlo generate A^ xpt events from the input 
energy spectrum, from which the maximum likelihood (hereafter 'ML') mass and cross- 
section are calculated. 

Finally we plot the distribution of ML masses and cross sections. Given real data 
from a single experiment a Bayesian analysis with priors on the WIMP parameters, 
possibly based on the results of other experiments, would be a reasonable approach. 
However, the question which we are trying to address ('For a given underlying WIMP 
mass, how well can the mass be measured?') is best answered by simply considering 
the distribution of ML masses and cross-sections. Since the approximations we make 
regarding the detector set-up and WIMP distribution and cross-section are optimistic, 
a real experiment would make less accurate determinations of the WIMP properties. 



3.4- Results 

In fig. [21 we plot the probability distribution of ML WIMP masses and cross-sections 
for input WIMP parameters m x = 100 GeV and a p = 10~ 7 pb. In this case both the 
input energy spectrum and the maximum likelihood analysis of the simulated events are 
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Figure 2. The distribution of the maximum likelihood WIMP masses, m x , and cross- 
sections, <7 p , for exposures of (top row, left to right and then bottom row left to right) 
£ = 3 x 10 2 , 3 x 10 3 , 3 x 10 4 and 3 x 10 5 kg day. For £ = 3 x 10 2 kg day we explicitly 
plot the results from all 10 4 Monte Carlo experiments. For the larger exposures we 
plot contours containing 68% and 95% of the probability distribution. In each panel 
the large cross denotes the input parameters: m x — 100 GeV, <r p = 10~ 7 pb. 



carried out assuming a Maxwellian speed distribution with v c = 220 km s" 1 . For each 
experiment the extended likelihood is maximized for WIMP parameters which produce 
an expected number of events equal to the actual number of events observed in that 
experiment: A(m x , <r p ) = A^ cxpt . This means that, for fixed exposure, the ML parameters 
are localized on curves corresponding to fixed N expt . For a given experiment the position 
of the ML parameters on the curve depends on the energies of the observed events. For 
£ = 3 x 10 2 kgday, A in = 7.8, which is sufficiently small that the stratification of ML 
parameters is clearly visible and we hence plot the actual pairs of m x — a p values. For 
the larger exposures the mean number of events expected is proportionately larger, the 
stratification is no longer visible, the ML values are better localized in the m x — a p plane 
and we instead plot contours containing 68% and 95% of the simulated experiments. 
We calculate the continuous probability distribution of m x and a p by smoothing the ML 
values from the 10 4 Monte Carlo simulations with a double gaussian kernel and summing 
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Figure 3. As fig. [2] for m x = 25 GcV. For £ = 3 x 10 2 kgday, A in = 4.2 and 
consequently there is a ~ 1% probability that an experiment would not detect any 
events. 

them. The deviations of the kernel are chosen, for each exposure and underlying mass, 
to produce relatively smooth contours without artificially inflating the spread in the ML 
values. We then integrate the probability density above a threshold value and vary this 
threshold to find the values which enclose 68% and 95% of the probability distribution, 
and plot contours corresponding to these threshold values. 

For the two smallest exposures there is a significant low probability density tail at 
large ML masses. This is due to the weak mass dependence of the characteristic energy 
Er for m x > rriA- For these exposures the majority of experiments have ML masses 
smaller than the input mass. The likelihood analysis is not biased, the expected value 
of the ML mass is equal to the input mass, however there are a significant number 
of experiments with ML masses substantially larger than the input mass. As the 
exposure is increased, Aj n becomes large and the fractional spread in the number of 
events observed in each experiment becomes small. This allows the interaction cross- 
section (which effectively acts as a normalization factor) to be accurately determined. 
The larger number of events also allows the energy dependence of the differential event 
rate, and hence the mass, to be determined more accurately. 
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Figure 4. As fig. [2] for m x = 50 GeV. 



With an exposure of £ = 3 x 10 2 kgday and an underlying WIMP mass of 
m x = 100 GeV it will be difficult to make quantitative statements about the WIMP 
mass, beyond excluding very large or very small masses. Similarly it will not be possible 
to measure the cross-section more accurately than to within an order of magnitude. For 
£ = 3 x 10 3 kgday the accuracy with which the mass can be determined improves, 
however there is still a significant tail of experiments finding ML masses substantially 
larger than the underlying value. More quantitatively 95% of experiments will find ML 
masses in the range 60 GeV < m x < 200 GeV. Increasing the exposure by an order of 
magnitude, the distribution of ML masses becomes more symmetric about the input 
mass and the spread in the ML masses is of order ±30 GeV. With a further increase 
in the exposure to £ = 3 x 10 5 kgday the mass could be measured to an accuracy of 
around 10 GeV (see the top panel of fig. [8] for a zoom in on the distribution in this case). 
The accuracy with which the cross-section can be measured increases with increasing 
exposure: A(loga p ) ~ ±0.2,0.05 and 0.02 for £ = 3 x 10 3 , 3 x 10 4 and 3 x 10 5 kgday 
respectively. The shape of the distribution of ML parameters reflects the shape of curves 
of constant A^ expt in the m x ~ 100 GeV region, with a weak positive correlation between 
the best-fit values of m x and <r p . 

In figs. [3]l6] we plot the results for input WIMP masses of m x = 25, 50, 250 and 
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Figure 5. As fig. [2] for 



250 GeV. For £ = 3 x 10 2 kgday, A in = 4.3 and 



consequently there is a ~ 1% probability that an experiment would not detect any 
events. For £ = 3 x 10 3 kgday there is a low probability density tail extending to 
large m x , and hence we can not accurately calculate a contour containing 95% of the 
probability density. The disconnected 'blobs' at large m x are also a consequence of 
this extended low probability density tail. 



500 GeV respectively. For m x = 25, 250 and 500 GeV the input mean number of events 
for E = 3 x 10 2 kgday is small enough that there is a significant probability that an 
experiment will see no events (making it impossible to determine the WIMP mass). 

For light WIMPs, m x < m^, the characteristic energy, varies significantly 
with WIMP mass. This allows the mass to be determined with large exposures more 
accurately than for m x = 100 GeV. For an underlying mass of m x = 25 GeV and an 
exposure of S — 2 x 10 3 kg day, due to the small expected number of events, it will be 
difficult to place meaningful constraints (beyond an upper limit on the mass) on the 
WIMP parameters. With larger exposures it will be possible to measure the WIMP mass 
and cross-section with increasing accuracy; for £ = 3 x 10 3 , 3 x 10 4 and 3 x 10 5 kg day the 
distribution of ML WIMP masses is symmetric and 95% of experiments lie within ±6, 4 
and 1 GeV of the input mass respectively. The accuracy with which a v can be measured 
improves roughly as for m x = 100 GeV. For m x ~ 50 GeV the fractional accuracy with 
which m x can be measured is similar; 95% of experiments lie within ±12, 7, 2 GeV of 
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Figure 6. As fig. [2] for m x = 500 GeV. For £ = 3 x 10 2 kgday, A in = 2.4 and 
there is a > 50% probability that an experiment would observe two or fewer events. 
Consequently we instead plot the 68% and 95% exclusion limits which would be found 
by an experiment observing N cxp = 2 events (the number of observed events with 
the largest probability). For this input mass there is a low probability density tail 
extending to large m x , even for £ = 3 x 10 5 kg day, and hence it is not possible to plot 
95% probability contours. 



the input mass for £ = 3 x 10 3 , 3 x 10 4 and 3 x 10 5 kgday respectively. The curves of 
constant N cxpt are, for m x = 50 GeV, roughly parallel to the m x axis, improving slightly 
the accuracy with which a p can be measured. 

For massive WIMPs, m x 3> m^, the weak dependence of on the WIMP mass 
means that there is a large spread in the distribution of ML masses even for large 
exposures. There is a small, but significant, probability that an experiment will happen 
to observe one or more unusually large energy events and hence find a large (> O(TeV)) 
ML mass. Because of the extremely weak dependence of the energy spectrum on the 
WIMP mass, the extended likelihood function also varies weakly with the mass in these 
cases (where the ML mass is of order a few TeV or greater, varying the WIMP mass 
by a factor of two from the ML value only changes the 4th significant figure of the log 
of the extended likelihood function). It is therefore not be possibly to determine the 
WIMP mass accurately in these cases. 
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For an underlying WIMP mass of m x = 250 GeV and an exposure of £ = 
3 x 10 2 kg day it will be difficult to place any meaningful constraints on the WIMP 
parameters. As the exposure is increased it will be possible to place a lower limit on 
the mass: for £ = 3 x 10 3 , 3 x 10 4 and 3 x 10 5 kg day, 95% of experiments have ML 
mass greater than 50, 125 and 200 GeV respectively. However it will not be possible to 
place an upper limit tighter than m x < 0(1 TeV) on the WIMP mass with a reasonable 
(> 68%) degree of confidence. The shape of the curves of constant iV ex p for m x ^> 
mean that it will also be difficult to constrain cr p . For £ = 3 x 10 4 and 3 x 10 5 kg, 
95% of experiments lie within A (log cr p ) ~ and ^; 4 of the input value of <r p (the 
asymmetric spread arises from the asymmetry in the distribution of ML masses). 

The situation is even worse for an underlying WIMP mass of m x = 500 GeV. For 
£ = 3 x 10 2 kg day there is a > 50% probability that an experiment will observe two or 
fewer events. In these circumstances rather than attempting to determine the WIMP 
mass and cross-section it would be more reasonable to instead determine the regions of 
WIMP mass-cross-section excluded. The observed number of events with the greatest 
probability is A^ exp = 2. In the relevant panel of Fig. [6] we therefore plot the exclusion 
limits from the 68% and 95% upper limits on the underlying mean number of events 
if this number of events were observed, A < 3.5 and 6.3 respectively. For £ = 3 x 10 4 
and 3 x 10 5 kgday, 95% of experiments have ML mass greater than 125 and 350 GeV 
respectively, however it will not be possible to place an upper limit on the WIMP mass 
at more than 68% confidence. Similarly it will only be possible to place a lower limit 
on Op. 

In fig. [7J we show the effect of the uncertainty in the value of the local circular 

speed, v c , on the determination of m x . For an input WIMP mass of m x = 100 GeV and 

an exposure of £ — 3 x 10 3 kgday we vary v c between 180 and 260kms _1 [27]. The 

likelihood analysis is however carried out assuming v c = 220 km s _1 . We see in fig. [TJthat 

there is a degeneracy between m x and v c . The kinetic energies of the incoming WIMPs 

depend on their mass and velocities. For larger (smaller) v c the incoming WIMPs have 

larger (smaller) mean kinetic energies than assumed, resulting in larger (smaller) ML 

mass values. This statement can be made more quantitative by differentiating the 

expression for the characteristic energy eq. (jHJ): 

Am v r , . Nn Af c , . 
x = _[i + ( mx/mA )]^ . (10) 

This gives, for an input WIMP mass of m x = 100 GeV a ~ 20 GeV shift in the WIMP 
mass from a 20kms _1 uncertainty in v c . However, as we see in fig. [7J the shape of the 
distribution of the ML parameters changes, in qualitatively the same way as when the 
underlying WIMP mass is changed. As v c is decreased (increased) the expected number 
of events increases (decreases) and hence the best-fit cross-sections typically decrease 
(increase). 

Finally in fig. [8] we examine the effect of the uncertainty in the detailed shape of the 
local velocity distribution. For an input WIMP mass of m x = 100 GeV and the largest 
exposure, £ = 3 x 10 5 kg day, we use as input the logarithmic ellipsoidal model (which is 



13 



-6.5 





100 200 300 400 500 



100 200 300 400 500 



o 





100 200 300 400 500 100 200 300 400 500 

m x (GeV) 



Figure 7. As fig. [2] but varying the input circular velocity v c = 180, 200, 240 and 
260 km s -1 with the exposure fixed at £ = 3 x 10 3 kgday. The maximum likelihood 
analysis is carried out assuming v c — 220 km s -1 . 



the simplest triaxial generalization of the standard isothermal sphere) [30] and two of the 
sets of parameter values previously considered in Ref. [21]: p = 0.9, q = 0.8, 7 = 0.07 
and p = 0.72, q = 0.7, 7 = 4.0 with the Earth located on the intermediate axis. These 
parameters correspond to axial ratios 1 : 0.78 : 0.48 and 1 : 0.45 : 0.38 (i.e. quite 
extreme triaxiality, especially in the second case) and, in both cases, velocity anisotropy 
(3 = 0.1. The local circular and escape speeds are kept fixed at the values used for 
the standard halo model. The likelihood analysis is carried out assuming the standard 
Maxwellian velocity distribution. Even for the second, quite extreme, model the shift 
in the typical ML mass from the underlying value is relatively small (less than 5 GeV). 
The small increase in the expected number of events for these halo models leads to a 
small downwards shift in the distribution of best-fit cross-sections. 
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Figure 8. As fig. [7] but for the standard halo model (top panel) and the two 
non-standard halo models described in the text with the exposure fixed at £ = 
3 x 10 5 kg day. The maximum likelihood analysis is carried out assuming the standard 
Maxwellian speed distribution. 



4. Validity of assumptions 

4-1. WIMP distribution 

Since the differential event rate is directly proportional to p x a x>^ an y uncertainty in the 
local WIMP density leads straightforwardly to an equivalent uncertainty in o~ p . We have 
fixed the local WIMP density to the 'standard value' of p x = 0.3 GeV cm" 3 . Refs. [311132] 
found, using various observations to constrain the parameters of a range of halo models, 
local densities in the range 0.2 — 0.8 GeV cm -3 , which would lead to a factor of a few 
uncertainty in the determination of the cross-section. 

We saw above that an uncertainty in the local circular velocity translates directly 
into an uncertainty in the measured WIMP mass. While the annual modulation signal 
depends sensitively on the WIMP velocity distribution [22], the time averaged differential 
event rates produced by smooth halo models are similar to that found under the standard 
assumption of a Maxwellian velocity distribution (e.g. [23J, [23] ) - Hence, as we saw above, 
the uncertainty in the speed distribution, for fixed v c , only leads to a small systematic 
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error in the determination of the WIMP mass (see fig. [H])- In other words, for smooth 
halo models the characteristic energy of the energy spectrum, and hence the WIMP 
mass determination, depends only weakly on the detailed kinetic energy distribution. 
These smooth halo models are derived by solving the collisionless Boltzmann equation 
which assumes that the dark matter distribution has reached a steady state. The 
assumption that the local dark matter distribution is dynamically mixed, and hence 
smooth, is, however, a strong assumption, which may not be valid [33| [31]. Helmi, 
White & Springel [35] found that the 'particles' in a simulated Milky- Way like halo 
in a ~ (kpc) 3 volume around the solar radius were relatively smoothly distributed. 
Direct detection experiments, however, probe the dark matter distribution on sub-milli- 
pc scales. The highest resolution simulations carried out to date have resolution of order 
(9(100 pc) and hence can not resolve the dark matter distribution on the relevant scales. 

In CDM cosmologies structure forms hierarchically and the local dark matter 
distribution will depend on the fate of the first, smallest, WIMP micro-halos to form. 
The mass of these microhalos is set by the WIMP microphysics in the early Universe [36J , 
specifically kinetic decoupling and free-streaming, and (depending on the WIMP's 
interaction properties) is expected to lie in the range 1O _12 M to 1O~ 4 M [37]. The 
dynamical evolution of these micro-halos is being studied |38j, however the detailed dark 
matter phase space distribution on sub-milli-pc scales is not yet known with any degree 
of certainty. If the local dark matter distribution consists of a small number of streams, 
with a-priori unknown velocities, then the energy spectrum would consist of a series of 
(sloping due to the energy dependence of the form factor) steps. The positions of these 
steps would depend on the WIMP mass and the (unknown) stream velocities, while the 
height of the steps would depend on the (unknown) stream density. In this case it would 
therefore not be possible to determine the WIMP mass from the energy spectrum. 

4-2. Negligible background 

From an experimental point of view, the most significant assumption is probably that 
of negligible background. Non-zero background could be incorporated (c.f. Ref. |39|) 
by simulating the recoil spectra produced by neutrons and including the background 
event rate (and additional parameters modeling the spectrum of the background events) 
in the maximum likelihood analysis. This would, however, require detailed modeling 
of the detector set-up and shielding. Adding additional parameters to the maximum 
likelihood analysis would clearly degrade the accuracy with which the WIMP mass and 
cross-section could be determined. The extent of the degradation would depend on the 
shape of the background energy spectrum (and how well it is known); the more similar 
it is to the WIMP spectrum the larger the errors will be. 

4-3. Other sources of systematic error 

Finite energy resolution and uncertainty in the form factor are other potential 
experimental sources of systematic error. The Helm form factor, with parameter values 
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as advocated by Lewin and Smith [7J, deviates by only of order 1% from that calculated 
using electron elastic scattering data [30]. We have checked that Gaussian energy 
resolution, with full width at half maximum of order 1 keV [H] , does not affect the 
WIMP parameters extracted from the energy spectrum 0. Both of these issues are 
therefore likely to be less important than non-zero background and/or fine-grained 
structure in the WIMP distribution. 

From a theoretical point of view the WIMP may have spin- dependent interactions El 
with the nucleon and/or different coupling to the proton and neutron (e.g. Ref. [3T1H2]). 
The measurement of the WIMP mass, in principle, in this case has been considered in 
Ref. [25]. 



5. Summary 

We have examined the accuracy with which it will be possible to determine the 
WIMP mass from the energy spectrum observed in a SuperCDMS-like direct detection 
experiment given optimistic assumptions about the WIMP distribution and detector 
properties. If the WIMP distribution is smooth, the differential event rate varies 
with energy, modulo the energy dependence of the detector form factor, as dR/dE oc 
exp (— E/E R ) where the characteristic energy, E R , depends on the WIMP mass, m x . 
For light WIMPs {m x <C where is the mass of the target nuclei) E R oc m 2 , 
while for heavy WIMPs (m x ^> ttia) -Er ~ const. Consequently for light WIMPs the 
energy spectrum is strongly dependent on the WIMP mass, allowing the mass to be 
measured fairly accurately. For heavy WIMPs the dependence on the WIMP mass is 
far weaker making it difficult to measure the mass. 

We have carried out Monte Carlo simulations of a SuperCDMS-like experiment 
composed of a Ge target with an energy threshold of 10 keV and zero background. We 
assumed, initially, that the local WIMP density is known and that the WIMP speed 
distribution is Maxwellian with local circular velocity v c = 220 km s _1 . For an optimistic 
interaction cross-section of cr p = 10 _7 pb, just below the current exclusion limits from 
the CDMS experiment [T2] , we considered a range of WIMP masses, 25 GeV < m x < 
500 GeV, and, efficiency weighted, exposures, 3 x 10 2 kgday < £ < 3 x 10 5 kgday. 
For £ — 3 x 10 2 kgday the expected number of events is small, the ML masses and 
cross-sections are stratified on curves of constant number of events and it would not be 
possible to obtain better than order of magnitude constraints on the WIMP parameters. 
For £ = 3 x 10 5 kg day and an input mass of m x = 100 GeV it would be possible, given 
the validity of the assumptions stated above, to measure the WIMP mass with accuracy 
±10 GeV and the fractional cross-section with accuracy A(loga p ) = ±0.02. The mass 
of lighter WIMPs could be measured more accurately, however for very light WIMPs, 

% Since the underlying differential event rate is, modulo the form factor, close to exponential, Gaussian 
smoothing only changes its shape for energies of order the resolution (which are below the energy 
threshold) . 

+ Natural Germanium contains 7.7% 73 Ge which is sensitive to spin dependent interactions [32] . 
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m x < 0(10 GeV), the number of events above the detector energy threshold would be 
too small to allow the mass to be measured accurately. For more massive WIMPs there 
is a significant tail of experiments with ML masses significantly larger than the input 
WIMP mass. For heavy WIMPs, m x > 0(500 GeV), even with £ = 3 x 10 5 kgday it 
will only be possibly to place lower limits on the WIMP mass and cross-section. 

We then examined the effect of varying the underlying WIMP speed distribution 
for an input WIMP mass of 100 GeV. A change of ±20kms _1 in the local circular 
speed, v c , leads to a shift in the distribution of best fit WIMP mass of roughly ±20 GeV 
(although an increase in v c leads, like an increase in the underlying WIMP mass, to a 
large tail of experiments with large best fit masses). Changing the shape of the WIMP 
velocity distribution, while keeping v c fixed, leads to only a small change in the input 
energy spectrum and hence the shift in the best fit WIMP masses is relatively small, 
< (9(5 GeV), even for quite extreme smooth halo models. There is, for smooth halo 
models, a factor of a few uncertainty in the local WIMP density [311 [32] which leads 
to a corresponding uncertainty in <j p . The assumption of a smooth WIMP distribution 
may well, however, not be valid. The local WIMP distribution, on sub milli-pc scales, 
may be composed of a (a priori unknown) number of discrete streams with unknown 
velocities. If this is the case it will not be possible to extract constraints on the WIMP 
mass from an observed signal. 

Finally we discussed the validity of the other assumptions made, namely negligible 
background, purely spin independent coupling, perfect energy resolution and known 
detector form factor. Of these negligible background is probably the most significant 
and could in principle be taken into account by simulating the recoil spectra produced by 
neutrons and including the background event rate and energy spectrum in the maximum 
likelihood analysis. 
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